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A fast acquisition algorithm for the Galileo suppressed carrier, subcarrier, and 
data symbol signals under low data rate, low signal-to-noise ratio (SNR) and 
high carrier phase-noise conditions has been developed. The algorithm employs 
a two-arm fast Fourier transform (FFT) method utilizing both the in-phase and 
quadrature-phase channels of the carrier. The use of both channels results in an 
improved SNR in the FFT acquisition, enabling the use of a shorter FFT period 
over which the carrier instability is expected to be less significant. The use of a 
two-arm FFT also enables subcarrier and symbol acquisition before carrier acquisi- 
tion. With the subcarrier and symbol loops locked first, the carrier can be acquired 
from an even shorter FFT period. Two- arm tracking loops are employed to lock the 
subcarrier and symbol loops with the carrier loop open. In addition, a new method 
is introduced for loop parameter modification to achieve the final (high) loop SNR 
in the shortest time possible. The fast acquisition algorithm is implemented in 
the Block V Receiver (BVR). This article describes the complete algorithm design , 
the extensive computer simulation work done for verification of the design and the 
analysis, implementation issues in the BVR, and the acquisition times of the al- 
gorithm. In the expected case of the Galileo spacecraft at Jupiter orbit insertion, 
Pd/Nq — H .6 dB-Hz, Rs ym = 16 symbols per sec, and the predicted acquisition 
time of the algorithm (to attain a 0.2-dB degradation from each loop to the output 
symbol SNR) is 38 sec. 


I. Introduction 

Due to the Galileo high-gain antenna failure, the spacecraft will be transmitting very low data-rate 
signals with very low signal-to-noise ratios (SNRs). For a significant portion of the mission, including 
the time of the Jupiter orbit insertion (JOI), the S-band (2.3-GHz) carrier signal will have high phase 
instability due to solar plasma process noise. Under these conditions, fast acquisition is extremely im- 
portant both to minimize data loss and because slower acquisition methods fail when the phase is not 
stable over the (long) acquisition time. An additional complication to the acquisition process is that, for 
the first time for the Deep Space Network, the received signal modulation type will be a fully suppressed 
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carrier with a fully suppressed subcarrier, which requires new algorithms. A fast acquisition algorithm 
for Galileo suppressed carrier, subcarrier, and data symbol signals under low data rate, low SNR, and 
high carrier phase-noise conditions has been developed. The algorithm is implemented in the Block V 
Receiver (BVR) [1] and is under consideration for implementation in the buffered telemetry demodulator 
(BTD). This article describes the complete algorithm design. Details of the various algorithms are given 
in the referenced internal memoranda and will be published in future issues of this report. 

A. The Approach 

The approach of the fast acquisition algorithm is to implement nonlinear detection algorithms to 
generate signals with sinusoidal components related to the subcarrier, symbol, and suppressed carrier 
signals, and then to estimate the signal frequencies and phases from the signals via fast Fourier transform 
(FFT) methods. This part of the fast acquisition algorithm will be referred to as the FFT acquisition. 
The corresponding tracking loops are then closed with initial frequencies and phases set to the FFT- 
estimated values. Each tracking loop is closed with a minimal initial loop SNR, and the final (high) 
loop SNR, necessary for negligible output-symbol SNR loss, is achieved by gradually narrowing the loop 
bandwidth (and transition window width if applicable) to final values resulting in the low SNR loss. (The 
transition window in the subcarrier and symbol loops is the window region about the data transition 
that is integrated [2], The window has a value in the range (0,1]. Windowing is implemented for noise 
reduction whenever square signals are tracked. Hence, windowing is used for subcarrier and symbol 
tracking, but not for carrier tracking.) This part will be referred to as the loop acquisition. Therefore, 
the fast acquisition algorithm can be described as a composition of two parts: the FFT acquisition and 
the loop acquisition. 

Multiple novel enhancement methods were incorporated into the fast acquisition algorithm design to 
acquire the subcarrier, symbol, and suppressed carrier in a much shorter acquisition time than if using the 
existing methods. In the FFT acquisition, a novel combination of the carrier and subcarrier in-phase (I) 
and quadrature-phase (Q) demodulated signals is utilized. This approach reduces the number of symbols 
required for the FFT at the required SNR by a factor of two. The new signals are a combination of all 
four II, IQ, QI, and QQ, where the II, IQ, QI, and QQ are the I and Q carrier and subcarrier demodulated 
signals' as shown in Fig. 1, instead of the traditional combination of only two signals. Using a complex 
FFT results in a single-tone, sign-correct frequency estimate. Therefore, further testing for the polarity of 
the frequency estimate is eliminated. In frequency estimation from the FFT, a sine interpolation method 
(Appendix) is used that results in refined estimates of the frequency and phase. 

In loop acquisition, a new approach to the entire loop acquisition process, from choosing the initial 
bandwidth and window values through narrowing these parameters at the fast rate, is introduced. Also, 
due to the design to initialize both the phase and frequency of the tracking loops with refined (FFT) 
estimates, the loops start out essentially in lock immediately after closure. An additional advantage of the 
loop starting in lock is that the bandwidth and window narrowing process can be activated immediately 
after loop closure. 

A new approach to suppressed carrier acquisition is employed when the carrier phase noise is high. 
In this case, the subcarrier and symbol are acquired first with the carrier loop open, which can be done 
because the subcarrier and symbol signals have lower phase noise due to the lower frequencies relative 
to the carrier. This approach of first acquiring the subcarrier and symbol is taken because, when the 
subcarrier tracking loop (SCL) and symbol synchronization loop (SSL) are locked, the SNR loss in the 
carrier acquisition FFT is significantly reduced, resulting in a shorter acquisition FFT time requirement 
during which the effects of high carrier phase noise are expected to be less significant. 

Two-arm SCL and SSL have been introduced for tracking the subcarrier and symbol with the carrier 
loop open. The two- arm loops (which use both the carrier I and Q arms) provide a 3-dB gain over the 
traditional one-arm loops in the low SNR case (to be reported on in a future article of this report). When 
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the carrier loop is closed, only additional noise is introduced by the use of the second arm; therefore, the 
conventional one-arm SCL and SSL are used. 



sin (£c T s m + % ) cos (o Jm+0 SC ) 


r (m) = \j~2P cos (w c / + 0 C + A d (0 sgn [sin ((o $c t + 0 SC )]) + n i n (0 1 f ^ 

Fig. 1. FFT-acquisition setup for the fast acquisition algorithm in the BVR. 

B. Outline 

In Section II, the fast acquisition algorithm process is described. In Section III, the FFT acquisition 
algorithm is described. The loop acquisition is described in Section IV. In Section V, the issues for 
implementing the algorithm in the Block V Receiver are discussed. Results of computer simulations are 
presented in Section VI. Finally, in Section VII, the acquisition times of the algorithm are presented. 

II. Algorithm Flowchart 

The flows of the algorithm are shown in Figs. 2 and 3 for two cases: when the carrier phase noise is 
high and when it is low, respectively. 

A. Algorithm Flowchart When Carrier Phase Noise Is High 

For a significant portion of its mission, the Galileo spacecraft is expected to have high carrier phase 
noise (e.g., the carrier phase and frequency can be moved significantly during the initial acquisition FFT 
time due to solar plasma). In this case, the subcarrier and symbol, which are relatively more stable due 
to lower frequencies relative to the carrier, are acquired first. By doing so, losses to the carrier acquisition 


85 













signal due to the subcarrier and symbol are reduced. One source of loss is from the fact that, before 
subcarrier acquisition, only the first harmonic of the subcarrier signal is captured for the FFT acquisition, 
compared to after subcarrier acquisition when all the harmonics are captured. Additional losses arise from 
the subcarrier and symbol phase errors. After reduction of the subcarrier and symbol losses, the SNR 
in the carrier acquisition is enhanced. As a result, a shorter acquisition time is required for the carrier 
acquisition FFT, during which the effects of the carrier dynamics are expected to be less significant. 


For this case, the flowchart of the acquisition algorithm is shown in Fig. 2. Step (1) will be to FFT- 
acquire the subcarrier. The subcarrier is acquired using N SCiaC q FFTs with N, Ci acq different, equally 
spaced symbol phase offsets. We use N ac<acq symbol phase offsets so that losses due to symbol phase 
error are minimized if the error in predicted symbol rate is small (Section III. A). In step (2), subcarrier 
acquisition is verified (see Section III.D). If the subcarrier was not acquired, a larger FFT size will be 
determined in step (9), and the acquisition process restarts from step (1)- 

After subcarrier FFT acquisition, a decision is made if the symbol rate is known accurately. For 
automating this decision, a method utilizing the relationship between the amplitudes of the N sc ,acq 
subcarrier acquisition FFTs is currently under development. If the symbol rate is known, symbol phase 
is estimated via interpolation of the subcarrier acquisition FFTs [step (4)]. If the symbol rate is not 
known, the symbol frequency and phase are estimated via an FFT approach [step (11)]. For the symbol 
acquisition FFT, the data stored from the subcarrier FFT acquisition are reused after refitting the stored 
data with the FFT-estimated subcarrier phase error. If a longer data set is required for symbol acquisition, 
the stored data set will be complemented with consecutively collected new data. By mainly using the 
stored data, little or no additional time is required for data collection [step (10)]. 


After estimation of the subcarrier and symbol frequency and phase via the FFT acquisition, the SCL 
and the SSL are initialized with the FFT-estimated frequencies and phases and closed. The bandwidth 
and window narrowing process is started immediately [steps (5) and (6)]. Note that the initial phase used 
to start the loop is the phase estimated for time t = t cUoop , the time at which the loop is closed. The 
subcarrier loop is closed with a square subcarrier reference; hence, all harmonics of the square wave are 
used after this point. Because the subcarrier and symbol are tracked while the carrier is open, two-arm 
SCL and SSL are used, which provide a 3-dB gain over the one-arm loops in the low SNR case. 

Next, the carrier is FFT acquired [step (7)]. New data can be collected for this step immediately after 
the subcarrier and symbol loops are closed, as these loops are then assumed to be in lock. Finally, when 
the carrier is FFT-acquired successfully, the carrier tracking loop will be initiated with the estimated 
frequency and phase and closed, and the bandwidth narrowing process will be started immediately [step 

( 8 )]- 


B. Algorithm Flowchart When the Carrier Phase Noise Is Low 

The flowchart of the acquisition algorithm when the carrier phase noise is relatively low, e.g., when the 
Sun-Earth-probe (SEP) angle is large, is shown in Fig. 3. In this case, the carrier frequency and phase 
are expected to be stable for the duration of the longer initial FFT time. (Recall that a long FFT time is 
required without the subcarrier and symbol lock, i.e., with the losses from using only the first harmonic 
of the subcarrier and from the subcarrier and symbol phase errors.) Then, the carrier and the subcarrier 
can be acquired simultaneously [step (1)]. (The acquisition signals for the carrier and the subcarrier are 
formed from analogous combinations of the same set of II, IQ, QI, and QQ and have the same detection 
statistics.) The rest of the flowchart is analogous to that for the case of a high carrier phase noise. Note 
that in step (9) both the carrier and subcarrier errors estimated from FFT acquisitions are refitted to the 
stored data from step (1). 
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Fig. 2. Algorithm flowchart when the carrier phase noise is high. 
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Fig. 3. Algorithm flowchart when the carrier phase noise is low. 


III. FFT Acquisition 

In this section, the FFT acquisition component of the fast acquisition algorithm is discussed. The 
setup of the receiver for the FFT acquisition is shown in Fig. 1. The input signal is modeled to be 

r(m) = \/‘iP cos (w c t + 9 C + A d(t) sgn [sin (w Jc t + 0 jc)]) T n in(0 \t=mT, (1) 
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where 


w c — 27 r/ c = carrier frequency 
O c = carrier phase 

A = 90 deg (modulation index for suppressed carrier) 
d(t) = data symbol of rate R sym symbols per sec 
u sc = 27t/ jc = subcarrier frequency 
9 SC = subcarrier phase 

«*n(0 = input noise, modeled as white Gaussian noise in analysis 
T 9 = sampling period of the signal 

m = sample index corresponding to a sampling period of T s sec 
The /(•) and Q(-) signals are defined as 

I(m) = r(m) cos (u c mT s + 0 C ) | lpf 
Q(m) ~ r(m) sin (u) c m7; + 0 C ) | LPF 

where l^pp indicates low-pass filtering, and 

w c = 27r/ c = predicted carrier frequency used in the carrier numerically controlled oscillator 

(NCO) 

9 C = predicted carrier phase used in the carrier NCO 
The II, IQ, QI, and QQ signals are defined as 


II(k) = I(m) sin (u>, c mT s + §, c ) | wF. y „ 

= AD{k) sin (&,«(£)) cos ( <j> e ,, c (k )) + N n (k) (2) 

IQ(k) = I(m) cos (u, c mT, + 9, c ) |/DF, ym 

= AD(k) sin ( <j> e ,c{k )) sin ( <j> e ,sc(k )) + Ni Q (k) (3) 

QI(k) = Q(m) sin {u }C mT, + 9 }C ) \iDF. yrn 

— AD(k) cos (<f> e c (k)) cos (<f> et , c (k)) + NQj(k) (4) 

QQ{k) = Q(m ) cos ( w fC mT , + 0 SC ) |/DF. ym 

= cos (<f>e,c(k)) sin (<p eiSC (k)) + N Q Q(k ) (5) 


where |/DF, ym indicates integration over the estimated symbol period £ is the sample index corre- 

sponding to the fcth overlapped integration period as shown for N SCi<lcq = 4 in Fig. 4(a), and 


89 



i nur: 


(a) 
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Fig. 4. Symbol acquisition: (a) overlapped integration intervals of II, IQ, Ql, and QQ for N sc acq = 4; 
and (b) relationship between <$> sym and the symbol integration interval ( / ) which results in the PSD 
with the largest peak ( P m ). 
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$e,c(&) — 
^e,5c(&) = 

A = 
D(k) = 


2nfsc = predicted subcarrier frequency used in the SCL NCO 

predicted subcarrier phase used in the SCL NCO 

carrier demodulation error 

subcarrier demodulation error 

—y/2P/ 7r 

the data symbol integrated over the estimated symbol period 


Note that to avoid significant degradation due to integration of the signals over the estimated symbol 
period, T sym (in the IDF sym in Fig. 1), the subcarrier and carrier frequency predict errors are required 
to be much less than the symbol rate. A rule of thumb is to require A/ carr _ or _, c < (R sym / 20) Hz, where 
A fcarr — or — sc is the maximum carrier or suppressed carrier frequency error. Note also that only the first 
harmonic of the square subcarrier is retained in the I’s and Q’s. 

The first letters in II, IQ, QI, and QQ indicate the in-phase (I) or quadrature-phase (Q) carrier 
demodulation, and the second letters indicate the I or Q subcarrier demodulation in the receiver. The 
definition of the I and Q demodulation (whether cosine or sine for I or Q) is arbitrary. A different 
definition of the I’s and Q s leads to a different combination of the I’s and Q’s in forming the acquisition 
signals. However, if from end to end the signal definition of the I’s and Q’s and the combinations of the 
I’s and the Q’s are consistent, and if the FFT estimated phase and frequency are applied to the loops 
with the correct polarization, the final FFT-acquisition result will be the same. 


The II(k ), IQ(k ), QI(k ), and QQ(k) are generated for N SCi acq overlapping, shifted integration periods 
°f T sym sec. The number of overlapped integrations per symbol, N sc>acq) is an arbitrary integer value 
greater than or equal to 4. In the Block V Receiver, the minimum required value of N sc>acq = 4 will be 
used. Hence, in the remainder of the article, N SC)Qcq = 4 is assumed. For this case, the overlapped symbol 
integration periods are shown in Fig. 4(a). 

The subcarrier, symbol, and suppressed carrier are FFT acquired (as shown in Fig. 1) using the four 
I’s and Q’s. The subcarrier, symbol, and suppressed carrier acquisition signals, X, c (), X 5J , m (*), and 
X corr (')> are formed from combinations of //(•), IQ ( ■), Q/(), and QQ ( ■). From each acquisition signal, 
a power spectral density (PSD) is formed as the magnitude square of the complex FFT of the signal, and 
the peak bin with PSD value above the detection threshold and bin location within the range of frequency 
search is detected (Section III .D) . From the peak bin and two bins about it, the frequency and phase 
(at t = 0) of the tone, /fft and <f>FFT are estimated using a sine interpolation method (Appendix). 
Combining FFT-estimated / fft and 4>fft an d predicted / and 0, refined estimates of the subcarrier, 
symbol, and suppressed carrier are obtained. 


In this section, for the I and Q demodulation setup of the BVR shown in Fig. 1, the acquisition signals 
are defined for the subcarrier, symbol, and suppressed carrier. For each FFT acquisition, the SNR in the 
acquisition FFT, SNRfft, defined as, 


SNR fft = 


tone power 
noise power 


= N data f ± N data ^-S L 


( 6 ) 


is presented. (Note that SNRfft can be defined as N^ ata times an equivalent symbol SNR, E*/N 0} 
which is the symbol SNR times a signal degradation factor, S L , in the acquisition signal.) Also for each 
case, the variances of the FFT-estimated frequency and phase are given. Finally, the determination of 
FFT parameters in the FFT acquisition is discussed. 
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A. Subcarrier FFT Acquisition 

Subcarrier acquisition is performed in step (1) of Figs. 2 and 3. At this point of the acquisition, the 
symbol frequency is assumed to be known, but the symbol phase is unknown. (Note that if A/Z jym , the 
error in the predicted symbol rate, is not negligible, the SNR of the acquisition FFT will be degraded, 
resulting in either the subcarrier being detected with a weaker SNR, or not being detected in the first 
FFT.) For the case in which A R sym is negligible, the integration of II, IQ, QI, and QQ over the estimated 
symbol period will have a constant symbol phase error, A (in cycles), during the FFT and contribute a 

loss in SNRfft' 


Dl = \[ 1 + Q-2A) 2 ] (7) 


where D\ is the expected symbol phase-error loss averaged over symbol probability. 

To minimize the loss due to D \ , the subcarrier acquisition signals will be formed for four ( N 9C>a cq ) 
symbol phase offsets, 1/4 cycle apart as A JC? o> A$ Cf i, A 5 Cj 2> <md X SCj 3. The SNRfft s °f ^5c,o> 

X sc>2 , and A% c>3 will be proportional to D and D\ 9> respectively, where A* (i = 0, 1, 2, 3) is 
the phase error between the actual input symbol epoch and the ith integration interval. The four symbol 
offset integration intervals and the A^s are shown in Figs. 4(a) and (b). Loss due to D 2 X is minimized in 
the subcarrier acqui sitio n by estimating the subcarrier and frequency from the PSD with the largest peak 
(minimum value of D x ), i.e., the largest peak out of the four PSDs is detected and the FFT containing 
the largest peak is used for FFT frequency and phase estimation, Jfft and <!>fft (at t = 0), respectively. 

The four subcarrier FFT acquisition signals, A r ac> *(n) (i = 0, 1, 2, 3), with the four equally spaced 
symbol offsets, are 


X sc ,i(n) = {H(i + 4n) + jIQ(i + 4n)) 2 +{QI(i + 4n) + jQQ{i + 4n)) 2 (8) 

where n = 0, 1, 2, • • • , N data - 1, and N da ta is the number of data samples in the acquisition FFT. Using 
Eqs. (2)-(5), these equations reduce to 


X ic ,i(n) = (AD(i + 4n)) V' 2 ^-‘( n) + n ic>i (n) 


(9) 


Each PSD of X sc ,i( n ) 0 = 0. 1. 2, 3) will have a peak tone at the bin closest to twice the error in the 
predicted frequency, i.e., the frequency estimated from the FFT (as in the Appendix), / F ft , is 

/ FFT ~ 2(/ sc - fsc) (I®) 


where f sc is the predicted subcarrier frequency. Similarly, the FFT-estimated phase, <{>fft (estimate of 
phase at t = 0, as in the Appendix), is doubled: 


<f>FFT ~ 2 (0 SC — 9 SC ) 


( 11 ) 


Therefore, the estimate of subcarrier frequency and phase are 
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( 12 ) 


fsc, FFT = fsc + 


Ifft 


<t>sc,FFT(t) = O sc + -{- 27T f SC} FFTt 


(13) 


This new method of combining all four II, IQ, QI, and QQ (using both I and Q channels of the carrier) 
in the FFT-acquisition signal, instead of the traditional use of only the carrier Q channel, results in a 3-dB 
improvement in the SNRfft • Hence, the number of data points necessary to attain a given SNR^ft 
is reduced by a factor of two. Furthermore, the frequency error estimate is sign correct, eliminating the 
need to further test for the polarity of the estimated frequency value. 

The SNR in each of the subcarrier acquisition FFTs is as follows: 1 


S N Rfft } { 


Ndata jy ^L,bin ' 




+ 2 {EJNo 


(14) 


where 


Ndata = number of symbols in the FFT 
^ = — ^—^(symbol SNR) 

= symbol integration loss in the ith FFT, Eq. (7) 

SL,bin — power loss due to the tone not falling exactly on an FFT bin 2 

_ _1 sin 2 (7rAfo n /a) ^ sin 2 (7rA^ n /a) 

N data sin 2 (7rA bin/(aN data )) ~ (7rA 6in /a) 2 

A ^ = fractional bin offset between the true tone bin and the peak FFT bin 
a = zero-padding factor in the FFT 
Sl,a/ = power loss due to the tone being away from dc 

_ sin 2 (7rA/ 5C T 5ym ) 

(irAf sc T sym ) 2 

Af sc = fsc fsc 


The estimated frequency and phase will have the following variances: 3 


1 The SNRs in the FFT acquisition of the subcarrier and carrier are derived in M. Aung, “Derivation of the SNR’s in the 
Subcarrier and Carrier FFT Acquisition in the Fast Acquisition Algorithm,” JPL Interoffice Memorandum 3338-94-054 
(internal document), Jet Propulsion Laboratory, Pasadena, California, June 6, 1994. 

2 The S Lfb i n , and FFT frequency interpolation algorithm (Appendix) are derived in S. A. Stephens, “An Analysis 

of FFT Tone Acquisition,” JPL Interoffice Memorandum 335.1-92-14 (internal document), Jet Propulsion Laboratory, 
Pasadena, California, May 14, 1992. 

3 The variances of the FFT-estimated frequencies and phases, the probability of FFT detection, and the detection threshold 
are discussed in S. Stephens and M. Aung, “FFT Estimation Variances,” JPL Interoffice Memorandum 335.1-94-DRAFT 
(internal document), Jet Propulsion Laboratory, Pasadena, California, August 1994. 
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FFT 


1 0.16 1 
4 T^ ft SNRfft 


Hz 2 


(15) 


t>,FFT 


(0 


= ( 2*) 2 * 


1 0.01345 2 

4 ~SN~Rfft + a J' FFT X 



rad 2 


(16) 


B. Symbol Acquisition 

After the subcarrier FFT acquisition, it is determined if the error in the predicted symbol rate error, 
A/?, ym , is negligible. If so, only the symbol phase is estimated from the subcarrier acquisition FFTs; 
otherwise, both the frequency and phase are estimated with a new FFT. 

1. Determination If Symbol Frequency Is Known Accurately. In step (3) of Figs. 2 and 3, 
a decision has to be made if the symbol rate is known accurately. If it is not known a priori, a possible 
method of making this decision is to compare the relative amplitudes of the N SCi acq subcarrier acquisition 
FFTs. The basic idea is that when the symbol frequency is known, the N SCiOC q FFTs will have amplitudes 
proportional to D\ f where A/s are phase errors l/N SC}acq cycles apart. On the other hand, when the 
symbol frequency is not known, the M SCt acq FFTs will have comparable amplitudes. (A method utilizing 
this principle is under development and will be reported on in a future article.) 

2. Symbol Phase Acquisition From Subcarrier Acquisition FFTs. When A R aym is negligible, 
symbol phase can be estimated from the four FFTs with A,' (i — 0, 1, 2, 3) as follows. In Fig. 4(b), the 
integration intervals of the input symbol stream, d(-), are shown for integration phase offsets of 0, 0.25, 
0.5, and 0.75 cycle (labeled interval index i =0, 1,2, and 3, respectively, and A* is the phase difference 
between the input symbol epoch and the t'th integration interval). Define m as the interval index of the 
integration interval that results in the PSD with the largest peak, P m . It can be seen that nr varies 
depending on the input symbol phase <f> sym (<^jym is defined with respect to the estimated symbol phase 
(j) sym = 0 at t = 0 at the beginning of the first symbol integration interval), as 


m = < 


r 3, 
2, 
1, 
U>, 


if 1/8 < <f> S ym <3/8 cycle 
if 3/8 < <p, y m < 5/ 8 cycle 
if 5/8 < <t> aym < 7/8 cycle 
otherwise 


Index m for various values of (p S ym is shown in Fig. 4(b), where the label P m is used to indicate the 
integration interval that results in the PSD with the largest peak. 

Hence, the symbol phase at t = 0 (with respect to the integration interval i = 0) can be estimated as 


^ym.FFrf 0) = ^1 “ — ) + A n 


(17) 


where A m is the estimate of A m , the phase error between the input symbol epoch and the integration 
interval 777. One derives A m by writing the peak amplitudes of the four PSDs (Pj, i = 0, 1, 2, 3) as being 
proportional to D\ , which is a quadratic function of A* [Eq. (7)]. Using the quadratic relationship plus 
the additional information that phase errors A,’s (i — 0, 1, 2, 3) are one-quarter of a cycle apart, A m is 
estimated as 4 * 


4 Details of symbol phase estimation and FFT phase estimation discussed in the Appendix are available in M. Aung and 

S. Stephens to E. W. Stone, “Fast Acquisition Algorithm Developed for the Block V Receiver Support of the Galileo Low 

Data- rate Mission,” JPL Interoffice Memorandum 3338-93-167 (internal document), Jet Propulsion Laboratory, Pasadena, 

California, August 25, 1993. 
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b + \/b 2 — 4ac 


(18) 


A 


m 


2 a 


where 


a — 4d 
b = -4 ( 1 


C= 1+1- 


N 

1 Y sc,acq 

2 


N., 


d = 


* jc,ac^ 
~ \/^t 




y/Pm + g + \JP m -g 

Pi = peak value of the /tli PSD 


g = floor 


N, 


sc,acq 


In summary, for this case, the estimated symbol parameters are 


f sym y FFT — Psym t pred 


<^ym,FFr(l) — 4>sym,FFT{ 0) + ^ Isym^FFT^ 


The variance of the estimated phase will be approximately 


(19) 

( 20 ) 


kFFT 


16 x SN Ftp ft 


rad" 


( 21 ) 


where SN Rfft is the SNR in the subcarrier acquisition FFT given in Eq. (14). 

3. Symbol FFT Acquisition When R S ym I s Unknown. When the symbol rate is not known, the 
symbol frequency and phase are estimated from one FFT of the acquisition signal. The acquisition signal 
is formed using the Fs and Q’s stored from the subcarrier acquisition, complemented with consecutively 
collected data if a larger number of symbols are required. Because most or all of the data are already 
available, little or no addition data collection time is required. Note from Eqs. (2)-(5) that the stored 
Fs and Q’s have subcarrier and carrier demodulation error, </> e>c and ^ jJC , respectively. Using estimates 
of these demodulation errors available from the subcarrier FFT acquisition (and carrier FFT acquisition 
if performed prior to symbol acquisition), the errors in the stored data are demodulated out before the 
data are used for the symbol acquisition signal. 

The overview of the algorithm is given here. 5 The approach is to form a signal X sym (n) (for n — 
0,1,2,.-., data ~ 1) proportional to D 2 (i + 4n) (i = 0, 1, 2, or 3) that can be shown to be periodic with 


5 Details of the algorithm and the SNR in the FFT are given in M. Aung to E. W. Stone, “FFT Acquisition of Symbol 
Frequency and Phase After Subcarrier Acquisition With the Carrier Loop Open in the Fast Acquisition Algorithm,” JPL 
Interoffice Memorandum 3338-94-075 (internal document), Jet Propulsion Laboratory, Pasadena, California, July 19, 1994, 
and in M. Aung to E. W. Stone, “FFT Acquisition of Symbol Frequency and Phase After Carrier and Subcarrier Acqui- 
sition in the Fast Acquisition Algorithm,” JPL Interoffice Memorandum 3338-94-076 (internal document), Jet Propulsion 
Laboratory, Pasadena, California, July 19, 1994. 
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a frequency equal to A R 3ym Hz [Fig. 4(a)]. Therefore, the FFT of D 2 {i + 4n) returns harmonic tones 
at lAR, ym (l = 0, 1, 2, • • ■), and the frequency and phase can be estimated (Appendix) from the first 
harmonic. To eliminate the sign ambiguity in the frequency estimate, a complex FFT is formed from the 
signal and its quadrature as 


— X $ym j{jl) "f~ j A 

For N sc>acq = 4, it can be shown [Fig. 4(a)] that 


( 22 ) 


Xsym t z( n ) — X S y m (n + 1 ) + jX S y m (jl) ( 2 ^) 

because the sequence X sym j(n) = X sym (n+ 1) is in quadrature with the sequence X syrn> Q(n ) = X 5ym (n), 
as they differ by 90 deg in phase, as shown in Fig. 4(a). 

A further modification is added to the algorithm where the acquisition signal is formed as A^y m (n) 
proportional to £)^(n). Then, note from Fig. 4(a) that the consecutive samples of Aj,y m (Ti) have over- 
lapping integration intervals. The outcome of overlapping is an improvement in the FFT SNR, and the 
estimated frequency and phase are 


f sym,F FT = fFFT ( 2 ^) 

<j>sym,FFT{i) = <f>FFT + 27T f 3 y m> FFTt “ ^ ( 25 ) 

where f FFF and <^fft( 0 are the frequency and phase estimates from the FFT of the acquisition signal 
X sym (k). 

4. Case A: Symbol FFT Acquisition, Retrofitting Only the Subcarrier. When the symbol 
FFT acquisition is performed after the subcarrier is acquired, with the carrier loop still open [Fig. 2, steps 
(10) and (11)], the estimated subcarrier frequency and phase errors are rotated out from each sample of 
the stored Fs and Q’s, and the acquisition signal is formed as 

X 3ym (n) = | Re{(Q7(u)+iQQ(n))e-^«-(")}+jRe{(77(n) + j/Q(n))e-^-^)} | 2 
= {$/(») COS (<^ e ,»c(«)) + QQ( n ) sin (<^e,jc(n))) 

+ {//(») COS (<t> e]SC (n)) + IQ{n) sin (<j> et9C (n)) } (26) 


= A 2 D 2 (n) cos 2 (<f>e iSC (n) - <j> e}SC {n)) 4- n 5y m ( FFT 
~ A 2 D 2 {n) -f n symiFF T 


where ^ e>5C ( ) is the estimate of the subcarrier demodulation error <j*e t 3 c(') in the stored II, IQ, QI, and 
QQ [Eqs. (2)-(5)], estimated from the subcarrier FFT acquisition. 
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5. Case B: Symbol FFT Acquisition, Retrofitting the Subcarrier and the Carrier. When 

symbol acquisition is performed after the subcarrier and the carrier are acquired [Fig. 3, steps (9) and 
(10)], the symbol frequency and phase are estimated in a manner analogous to case A, with the following 
acquisition signal. The estimated subcarrier And carrier errors are rotated out from each sample of the 
stored I’s and Q’s, and the acquisition signal is 


= (Re [e-^'‘ (n) [Re {(Q/(n) + jQQ{n))e~>+‘"^} + jRe {(//(n) + 


= ( QI(n) cos + QQ(n) sin (* e, 3C (n))j cos (<^e,c(»)) 


+ 


|7/(n) cos (<£ e , «(»))+ IQ(n) sin (<^ eiJC (ra))j sin (0 e ,c(rc))) 


(27) 


= A 1 2 D(n) 2 cos 2 (4>e, s c(n) ~ 4>e,»c(n)) COS 2 (^ )C - <£ e ,c) + n»jm,FfT 
Ki A 2 D(n) 2 + n, ymt FFT 


where<£ ejc (-) and <fi e> c(‘) are the estimates of the subcarrier and carrier demodulation errors in the stored 
data, estimated from the subcarrier and carrier FFT acquisitions. 

For cases A and B, the estimated frequencies and phases will have the following variances: 


2 _ 0.16 1 
LFFT ~ T 2 ft SNRfft 


a 


(28) 


— ( 2?r ) 2 a f- FFT x ^p") ^ ra< ^ 2 (29) 

C. Carrier FFT Acquisition 

The suppressed carrier is acquired after the sub carrier and symbol acquisition when the carrier phase 
noise is high and simultaneously with the sub carrier when the carrier phase noise is low. The carrier 
acquisition is described for the two cases. 

1. Case A: Carrier Acquisition After the Subcarrier and Symbol Have Been Acquired. 

If the carrier is acquired after the subcarrier and symbol loops are in lock [Fig. 2, step (7)], the carrier 
frequency and phase are estimated from one FFT. Note that at this point of the acquisition process, the 
I and Q subcarrier demodulation is performed with square-wave references. The acquisition signal is 


Acarr(^) — (Q/(4n)+j//(4n)) 2 


2 


+ ™carr } acqijl) 


(30) 


The SNR in the FFT will be 


SNRfft = Nd aia —SL,bin — 
No 8 


02 


[•Sl.A/ + 2(EJN 0 ) 


(31) 
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2. Case B: Simultaneous Acquisition of the Carrier and the Subcarrier. For simultaneous 
subcarrier and carrier acquisition [Fig. 3, step (1)], the four carrier FFT acquisition signals are analogous 
to the subcarrier acquisition signal: 

X carrti (n) = ( QI(i + 4n) + jll{i + 4n)) 2 + (QQ(i + 4n) + jIQ(i + 4n)) 2 

= (AD(i + 4n) V 2 ^' c(n) ) + n carr: acq(n) (32) 

where i = 0, 1, 2, 3, and n = 0, 1, • • •, N data - 1. The SNR of the tone in the FFT and the variances of the 
estimated frequency and phase are the same as those for subcarrier acquisition, given in Eqs. (14)— ( 16) . 

Analogous to the subcarrier FFT acquisition, the estimated carrier frequency and phase for both cases 
A and B are 


fc,FFT 

i , Ifft 

- /c + 2 

(33) 

<f>c,FFT{t) 

= 0 C + 2 + 27T fc,FFTt 

(34) 


where f FFT and 4> FF j are the estimates of frequency and phase (at t — 0) of the peak tone in the FFT 
(Appendix), and f c and <f) c are the predicted carrier frequency and phase. The estimated frequency and 
phase will have the same variance as in Eqs. (15) and (16). 

D. Determination of the FFT Parameters 

In this section, the FFT parameters, namely Ndata, detection threshold, and the zero-padding factor, 
are discussed. 

1. Determination of Ndata- The number of data points to be used in the FFT, Ndata , is determined 
from choosing a desired confidence of detection, (7, and solving for Ndata • The confidence of detection, 
which is the probability of detecting the correct signal bin in the FFT in the presence of noise, is 



0.5e 


-0.5 (E:/N 0 )N daia 


2( A/, ea rch 


(35) 


where El/N 0 is the equivalent symbol SNR in the acquisition FFT, defined in Eq. (6). The Af sea rch 
determines the frequency search region, ±A /search Hz about the predicted frequency, i.e., the peak bin 
in the PSD is detected as the bin with the largest PSD value among the bins corresponding to the search 
frequency range [/ — Af sea rch, f + A/ search] Hz. It should be noted that Eq. (35) is only an approximation 
of the confidence of detection. 

2. Detection Threshold. In the FFT acquisition, FFT detection is declared if the peak bin in the 
search frequency range [/ — A f sea rch,f + search] Hz exceeds a lower threshold, Pthrcsh * The lower 

threshold is set as 

Pthresh = (P, ( 36 ) 


where 
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P s = the expected signal power 

i= (l — Erf~ 1 [2C — l]) 2 

C = confidence of detection 


Erf- l (-) 


inverse error function 


3. The Zero-Padding Factor. The nominal value of the zero-padding factor, a, will be 4. However, 
Ndata is chosen according to Eq. (35), which may not be a 2 n value. In general, the FFT size, Nfft , 
will be determined as 


Nfft — 4 x N neare3 t 


where N nearest is the nearest 2 n value greater than or equal to Ndata - Data points Ndata plus ( Nfft 
— Ndata) zeros will be used to form the Nfft point FFT. Therefore, the actual zero-padding factor, 
which should be the value used in all equations requiring a (e.g., frequency interpolation), will be 


Nf ft 
Ndata 


> 4 


(37) 


IV. Loop Acquisition 

This section discusses loop acquisition. First, initialization of the tracking loop frequency and phase 
is described. The use of two-arm versus one-arm loops is also briefly discussed, and the majority of the 
section describes the process of initializing and gradually narrowing the loop bandwidths (and windows 
if applicable) until the final loop SNR is achieved. 

A. Initializing the NCOs of the Tracking Loops 

After the FFT acquisition, the (subcarrier, symbol, and carrier if applicable) tracking loops are closed 
with the FFT estimated parameters. In each loop (indicating subcarrier, symbol, or carrier as Y = 
“sc,” “sym,” or “carr” in the following), the initial frequency in the numerically controlled oscillator 
(NCO) is set to fNCOjnit' 


f NCO t init = fv,FFT 

and the initial phase of the NCO is set to <j>NCO,in%F 


(38) 


<f>NCO ,init = <f>Y,FFT(tclJoop) 


(39) 


where t = tdjoop is the time at which the loop is closed with respect to / = 0 at the beginning of the first 
symbol integration, shown in Fig. 4(a). 

B. Two-Arm and One-Arm Tracking Loops 

New two-arm tracking loops have been introduced for subcarrier and symbol tracking while the carrier 
loop is open. The two-arm loop forms the loop phase error estimate (f> e> loop from two phase detector 
outputs, one output using the carrier I input and the other output using the carrier Q input, as shown 


99 



ill! till II l!ll I 


in Fig. 5. When tracking the subcarrier and/or symbol loop while the carrier is open-loop demodulated, 
and when the symbol SNR is low, using two arms provides a 3-dB gain in the loop SNR as the signals in 
both channels are used. When the carrier is in lock, no additional information and only additional noise 
is gained by using the second channel and, therefore, a conventional one-arm loop is used. 


(a) 



Fig. 5. Subcarrier or symbol tracking loops: (a) one arm and (b) two arm. 


C. Initial Bandwidths and Windows 

When FFT acquisition is complete, each estimate of signal phase and frequency will have associated 
with it a variance, FFF and cr'j respectively (given in the last section). The initial tracking loop 
bandwidth (and transition window if applicable) is chosen by matching the loop phase variance, 
with the FFT-estimated phase variance, v \ iFFT . This is analogous to a constant-bandwidth tracking 
loop handing over to the acquisition loops, which will start at the bandwidth and window and begin 
narrowing. 

The tracking variance a % for each type of loop can be written in the form 


(7 


2 


(2tt ) 2 c^Blw rad 2 , for loops with transition windows (SCL and SSL) 
(27 t) 2 c^Bl rad 2 , otherwise (suppressed carrier loop) 
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where B L is the loop noise bandwidth, w is the window width size, and is a constant that depends 
on the signal SNR and the tracking loop. For loops with transition windows, the window is specified in 
terms of the tracking variance with the parameter c w : 


Cw O (k 

w = 

_ f 4 c f wl for subcarrier loops 
\ 2 <4, for symbol loops 

The parameterization of c w in terms of <4 is f° r physical significance: <4 gives the tracking error that 
reaches the edge of the transition window, expressed in tracking standard deviations. For example, with 
<4 = 4, the (unlikely) tracking error of 4cr^ is necessary to cause the signal transition to fall outside of 
the phase extractor window. 

We can now solve for the initial bandwidth and window width: 


w(0) 


CwO^^FT 

2w 


Bl( 0 ) 


2 

a <j>,FFT 

(2tt) 2 cI w (Q) 


For a loop with no transition window, the initial bandwidth is 


Bl( 0) = 


2 

^<p,FFT 

(2*) 2 4 


The windowing allows the initial bandwidth to be w times larger than with no windowing for a given 
loop SNR. 

D. Bandwidth and Window Narrowing 

In the continuous bandwidth and window narrowing methods used, 6 the loop time constant r = 
1/(2B l ) is linearly increasing in time with rate dr/dt. Bandwidth as a function of time then has the 
form 


BlU) = 


1 


,2B L (0) 



(40) 


In tracking loops with transition windows, window as a function of time is given by 


w(i) = u' (0) 


Bdt) 
Bl{ 0) 


6 The bandwidth and window narrowing method is introduced and analyzed in S. Stephens, “Bandwidth and Window 
Modification in Phase- locked Loops,” JPL Interoffice Memorandum 335.9-94-010 (internal document), Jet Propulsion 
Laboratory, Pasadena, California, June 1994. 
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From simulation 7 and analysis, values of dr/dt and c' w that minimize the loop phase variance after a fixed 
time span are specified for loops 8 of order N and damping parameter t / 2 and are given in Table 1. 


Table 1. Optimum values of dr/dt and d w for 
loops of order N and damping parameter r; 2 . 


Loop 

dr/dt 

4. 

Subcarrier loop, N = 2,tj 2 = — 1 

0.181 

4.0 

Symbol loop, N = 1 

1.0 

4.0 

Symbol loop, N = 2, r? 2 = —1 

0.181 

4.0 

Carrier loop, N — 2,?j 2 = — 1 

0.383 



From solving Eq. (40), the time, Tji n , to reach from Bl(0) is 

T = _L f-L l —) L_ 

/,n 2f \B L j in B L (0)J ~ 2 %B LJil 


Note that for dr/dt — 0.181, 


(41) 


Tji 


2.8 

BlJiti 


sec 


For digital loops, bandwidth and window width values can be changed no faster than the update 
rate. Processor limitations may impose an even tighter constraint: for the Block V Receiver, the loop 
parameters may only be modified once per second. If the continuous equations above are used to generate 
discrete values of bandwidths and windows at rates ( f u ) much larger than the loop bandwidth, then the 
analysis should still hold; this is confirmed by simulation in Section VI. However, if f u is comparable to 
or larger than the loop bandwidth, the analysis in Footnote 6 specifies a maximum allowable change in 
the value per update time of 1 j f u sec. In this case, the bandwidths and windows are narrowed in the 
limited step sizes allowed per update time. For Galileo fast-acquisition in the BVR (/ u — 1 Hz), where 
the bandwidths are less than f u = 1 Hz for all loops, we can always use the continuous approximation. 

The bandwidth function versus time is illustrated in Fig. 6(a). The dark line shows the actual band- 
width values applied to the tracking loop. The steps reflect the update rate f u = 1 Hz, at which the 
bandwidth is modified. The following sections discuss BL >m in and Biji n . 

The window function versus time is illustrated in Fig. 6(b). The dark line shows the actual window 
width values applied to the tracking loop. The limitation of the window value (in the BVR) to a 2 n 
value (discussed in Section V) results in quantized window values, as in the figure. A discussion of wji n 
and Wj in follows. 


7 The implementation and simulation of Footnote 6 is discussed in M. Aung, S. Stephens, and C. Buu to J. Berner, 
“Block V Receiver Implementation and Simulation of Bandwidth and Window Narrowing Methods in [l],” JPL Interoffice 
Memorandum 3338-94-042 (internal document), Jet Propulsion Laboratory, Pasadena, California, June 3, 1994. 

8 S. A. Stephens and J. B. Thomas, “First-Principles Analysis of Digital Phase-Locked Loops with Controlled-Roots, 
Transient-Response-Specific Parameterization,” JPL Interoffice Memorandum 335.1-91-032 (internal document), Jet 
Propulsion Laboratory, Pasadena, California, November 22, 1991, 
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(a) B L {f) 



(b) co (0 



Fig. 6. Bandwidth and window narrowing: (a) loop bandwidth versus time fand (b) co (t ) versus time f. 


E. Final Bandwidths and Windows 

The final bandwidth, is determined from a specified maximum output symbol SNR loss; the 

bandwidth narrowing routine stops narrowing when it reaches the final required bandwidth at time 
t — Tji n . However, window quantization will cause the final quantized window value u>j in to be larger 
than wf in , specified by the continuous equations above. This in turn yields a lower final loop SNR 
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than desired. For compensation, the window or bandwidth must be narrowed an additional amount. 
During narrowing, the tracking variance at any given time can be twice as large as that implied by the 
bandwidth alone, due to transients in the narrowing process. Thus, the optimum windowing factor d w 
during narrowing will be larger than the factor used once narrowing transients have died away. The 
method implemented for reaching the final design loop SNR makes use of this smaller c' w when the 
narrowing process is completed; the window is narrowed by an additional factor of two at the end of the 
narrowing process, which is when the final bandwidth is reached,, as shown in Fig. 6(b), i.e., at t Tjxm 
the window is set to w'j in /2. 

F. Special Considerations 

For the symbol loop, the initial bandwidth specified above may cause the initial loop SNR of the SSL 
to be somewhat low (e.g., when S N Rfft = 12 dB and symbol phase is estimated from a set of subcarrier 
acquisition FFTs (Section III.B.2), the initial SSL SNR is about 8 dB). For this case, the bandwidth is 
limited to Bt.min, which corresponds to a loop SNR of 10 dB; the window width calculation uses the 
continuous bandwidth function before it is limited, so that the window is kept larger for transients caused 
by the narrow bandwidth [Fig. 6(a)]. 


V. Implementation in the Block V Receiver 

The fast acquisition algorithm has been implemented in the BVR; this section discusses the issues 
regarding its implementation. The BVR high-speed digital signal-processing hardware consists of five 
application-specific integrated circuits (ASICs) that are used to demodulate the carrier and up to two 
subcarriers and two symbol streams. There are three unique ASIC designs: carrier processing, subcarrier 
processing, and symbol processing. The two telemetry processing channels (subcarrier and symbol loops) 
are provided by independent sets of the subcarrier and symbol ASICs, each set receiving the in-phase (I) 
and quadrature-phase (Q) carrier demodulation outputs (the subcarrier ASIC can internally reverse the 
I and Q inputs). All tracking loop NCOs are frequency and phase updatable. The software processing of 
each telemetry channel is also independent of the other. 

Implementation of the fast acquisition algorithm involved changing the software tasking structure of 
the BVR subcarrier and symbol tracking loops and modifying the algorithm for the software implementa- 
tion. The carrier, subcarrier, and symbol tracking loops are implemented as independent software tasks. 
As such, the loop acquisitions are independent processes; it is possible to acquire the subcarrier and sym- 
bol loops sequentially or simultaneously. The fast acquisition algorithm requires that the subcarrier and 
symbol loops be controlled together. To accomplish this, a separate task was created for the fast acqui- 
sition process. The subcarrier and symbol loops have two acquisition tasks each, the normal acquisition, 
which is independent for each loop, and the fast acquisition, which controls both loops. 


Upon receipt of the acquisition command, the task zeroes out the subcarrier and symbol NCO phases 
and sets the subcarrier NCO frequency to the predicted subcarrier frequency and the symbol NCO 
frequency to the predicted symbol rate times V 3C , ac? . The subcarrier reference signal type is set to sine 
wave (as opposed to the normal square wave). Because there are no accumulators to sum the symbol 
integration outputs, the symbols must be read into software for processing. This means that the symbol 
NCO frequency used must be less than the symbol hardware interrupt rate (1 or 2 kHz). This places a 
limit on the product of V jc , ac , and the symbol rate. For this part of the acquisition process, only one set 
of telemetry ASICs is used, since all four outputs (II, IQ, QI, and QQ) are available. 

The FFT size is determined from the input predicted signal parameters. The data are collected, and 
the complex subcarrier FFT inputs for the N SCiacq FFTs are generated. A Sky Computer Skybolt array 
processor board is used to compute the FFTs, which can be as large as 2 17 points. After the (subcarrier, 
symbol, and carrier) phase and frequency are estimated from the FFT acquisition, the tracking loop 
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NCOs, which are frequency and phase updatable, are initialized with the FFT-estimated values, and the 
loops are closed. The subcarrier reference signal type is changed to the square wave. 

If the carrier loop has not yet been acquired, the subcarrier and symbol loops must be operated in the 
two-arm mode. This uses both sets of telemetry ASICs to generate the windowed subcarrier and symbol 
tracking loop inputs (the second set of ASICs has the I and Q inputs reversed, providing the needed 
signals to give windowed versions of the II, IQ, QI, and QQ signals). If the carrier is locked, only one 
telemetry channel is required. 

The implementation of the transition window-narrowing and bandwidth-narrowing algorithms also 
required some trade-offs, due both to hardware and software concerns. The subcarrier and symbol 
windows are implemented in the ASICs and are quantized to a factor of two. In other words, the 
hardware window widths are step selectable in increments of 2 _n , where n is an integer in the range of 0 
to 15. So, the algorithm could not continuously vary the window width; it had to quantize the windows 
to the hardware capabilities. 

The update rates of the narrowing algorithms (both window and bandwidth) were also quantized in 
time. All tracking loops have two update rates: the hardware interrupt rate, which is the rate that the 
software reads in the data off of the ASICs and processes the tracking loop, and the once-per-second pulse, 
which is the rate at which the performance data are reported and the loop configuration parameters are 
updated. It was determined that the narrowing algorithms could be updated at the lower rate (once per 
second), which allowed a simpler implementation that did not require lowering the interrupt rate (the 
higher the interrupt rate, the less time there is to process the software instructions required to close the 
loop). 


VI. Computer Simulation 

Extensive computer simulations were performed to verify the design and analysis of the fast acquisition 
method. 9 The simulations examine the following aspects of the fast acquisition algorithm: 

(1) Subcarrier FFT acquisition 

(2) Symbol FFT acquisition 

(3) Carrier FFT acquisition 

(4) Subcarrjer, symbol, and carrier loop acquisition 

In order to verify case 1 [step (1) of Figs. 2 and 3], the subcarrier FFT acquisition is simulated 
in the presence of demodulation errors, and the resulting frequency and phase estimates are compared 
to the input parameters. The theoretical performance of the algorithm is also evaluated in terms of 
the probability of FFT detection and statistics on the frequency and phase estimates. Specifically, the 
simulations verify agreement between measured and predicted values for the following: 


(la) Subcarrier tone SNR when the symbol rate error is negligible and the subcarrier is acquired 
using four initial symbol phase offsets 

(lb) Subcarrier tone SNR when the symbol rate error is significant 


9 References for the simulation results are summarized in M. Aung to E. W. Stone, “Computer Simulation Work Sup- 
porting the Design and Analysis of the Block V Receiver Fast Acquisition Algorithm,” JPL Interoffice Memorandum 
3338-94-053 (internal document), Jet Propulsion Laboratory, Pasadena, California, June 1, 1994, which gives an overview 
of the simulation work and lists all interoffice memoranda reporting the simulation results. 
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(lc) FFT detection confidence 

(l d) Subcarrier frequency and phase estimates 

(le) Subcarrier frequency and phase estimate variances 

Verification of (le) is important since the SCL is closed with the initial bandwidth, and the window values 
are calculated as functions of the predicted standard deviation of the FFT-estimated phase. In addition, 
(la) and (Id) were verified for actual Galileo downlink data recorded at Goldstone. Good agreement to 
within one or two standard deviations of the measured value was seen for (la), (lb), (Id), and (le). For 
(lc), an approximate predict was used, resulting in some divergence with the simulated value. 

In symbol FFT acquisition, when the error in the predicted symbol frequency is negligible, the symbol 
phase can be estimated from the subcarrier acquisition by interpolating the N SCi acq PSDs of the subcarrier 
FFT acquisition If the symbol phase detection fails due to a symbol rate error, the symbol frequency 
and phase are acquired via FFT methods. Hence, for steps (4), (10), and (11) of Fig. 2 and steps (4), 
(9), and (10) of Fig. 3, the simulations verify agreement between measured and predicted values for the 

following: 


(2a) Symbol phase when the symbol rate error is negligible and the phase is estimated from the four 
subcarrier acquisition FFTs 

(2b) Symbol phase estimate variances for (2a) 

(2c) Symbol frequency and phase and the variances of these estimates when the symbol rate is 
unknown 

(2d) Symbol error tone SNR when the symbol rate is unknown 

Verification of variance of the FFT estimate is important since the SSL is closed with the initial bandwidth 
and window values, which are calculated as functions of the predicted standard deviation of the FFT- 
estimated phase. For (2a) through (2d), good agreement to within one or two standard deviations of the 
measured value was seen. Agreement varied slightly for (2b), but the uncertainty lies in the computation 
of the predicted standard deviation (an approximation of the true standard deviation was used). 

Simulations of the subcarrier and carrier simultaneous FFT acquisition [step (1), Fig. 3] confirm that 
carrier FFT acquisition is analogous to subcarrier FFT acquisition in terms of tone SNRs. Therefore, 
results obtained for case (1) can be applied to case (3). 

The loop acquisition in case (4) was simulated separately by initializing the loop NCO phase and 
frequency with errors equal to one standard deviation of the FFT estimates. The loops are then closed 
with initial bandwidths and window widths as described in Section IV.C. The band widths are narrowed 
as described in Section IV. D until the desired symbol SNR degradation is achieved. The window widths 
are narrowed in quantized increments of 2“" to reflect hardware implementation constraints in the BVR. 
In addition, the simulation update rate of the narrowing algorithm is set to 1 pulse per sec in accordance 
with the BVR loop configuration parameter update rate. The simulations measure acquisition times for 
the following cases: 


(4a) The subcarrier, symbol, and carrier loops are closed simultaneously using the one-arm approach 
[steps (5) through (7) in Fig. 3]. 
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(4b) The two-arm SCL and SSL are closed, followed by carrier FFT acquisition first and then the 
carrier loop acquisition [steps (5) through (8) in Fig. 2]. This is simulated by closing all three 
loops simultaneously, assuming that the carrier FFT acquisition time is negligible compared to 
the SCL and SSL loop acquisition times. 

(4c) The two-arm SCL and SSL are closed while the carrier loop is left open [steps (5) and (6) in 
Fig. 2], 

Good agreement is observed between measured and predicted acquisition times as discussed in Section VII. 


VII. Acquisition Times 

The acquisition time of the fast acquisition algorithm is defined as the length of time required until 
all the loops (subcarrier, symbol, and carrier if applicable) have attained loop SNRs high enough to 
contribute only a tolerable amount of degradation to the output symbol SNR. For the high carrier phase- 
noise case (Section II.A) when the carrier is acquired after the subcarrier and symbol acquisition, the 
acquisition time is 


T acq = MAX{T, c>acq 


L sym,acq 


}+T c 


carr,acq 


sec 


(42) 


where Ty jQcq ( Y = sc , sym , or carr , indicating subcarrier, symbol or carrier) is the acquisition time of 
the indicated signal. For the low carrier phase-noise case (Section II. B) when the subcarrier, symbol, and 
carrier are acquired simultaneously, the acquisition time is 


T 

acq 


M AX\T sc>acq 


i 'Rsym,acq > 


T 

-*■ carr, acq 


} 


sec 


(43) 


In the following, individual acquisition times for the subcarrier, symbol, and suppressed carrier 
{T sc ,acq> T s y mt ac q > and T carrtacqj respectively) are presented. For each case, Ty acq is defined as 


Ty t acq = TV, FFT -f Tyjoop 


(44) 


where Ty )FFT is the FFT acquisition time and T Yf i oop is the loop acquisition time. 
The FFT acquisition time, TV,fft , is 


7y, FFT = N daia T $ym = T$ y m (45) 

where SNR Y FF t is the SNR in the FFT-acquisition of the signal under consideration. 

The loop acquisition time for each loop, Tyj oop , is defined as the length of time for narrowing the 
bandwidth and window width until the design loop SNR, SNR L j in , is achieved. The SNR LJin is 
determined by the maximum allowed degradation, Deg(<f>) y to the output symbol SNR due to the loop 
phase jitter, <j> . The relationships between Deg(<£) a nd <j> are listed in Table 2 for the SSL, SCL, and the 
suppressed carrier loop. The expected degradation, Deg(<j>) t calculated as the statistical expectation of 
Deg(4>) assuming a Gaussian distribution for <f>, is also listed in Table 2 where SNR L = l/cr*. 
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Table 2. Degradation, Deg(<f>) t to the output s ymbol S NR due to loop phase 
error <f>, and expected values of Deg(4>) and Deg(<j>) for the SCL, SSL, and 
carrier loops. 


Loop 

Deg(4>) 

Deg(4>) 

SCL 

(-<£)’ 

1 1 4 

1 y 2SNR l tt 3 + snr l * 2 

SSL 

i(i + (i- 2 M) 2 ) 

i 1 2 

1 4 y SNR l { 2») 3 + SNR l { 2tt)2 

Suppressed carrier loop 

COS 2 (4>) 

1 (l + e~ si ? h l ^ 


The SNR L jin is calculated as the minimum SNR L required for a given expected degradation. In 
Table 3, SNRLjin is calculated for expected degradation values of { 0 . 1 , 0 . 2 , 0 . 5 , 1 } dB, i.e., for final 
degradation of { 0 . 1 , 0 . 2 , 0 . 5 , 1 } dB to the output symbol SNR due to each loop, the loop bandwidth 
and window width are to be narrowed until the loop SNR of SN RLjin Table 3 is achieved. When 

designing for low final loop SNRs, the SN RljIu derived should be padded conservatively, e.g., by 1 dB, 
to compensate for nonlinear effects not modeled in the equations used for loop SNRs, 

Table 3. Minimum required loop SNRs, SNRi ffn s, as a 
function of Deg(<f>). 


Deg(4>)y 

dB 

SCL SNR L jin , 
dB 

SSL SNR L j in , 
dB 

Suppressed carrier 
loop SNRLjin i 
db 

0.1 

33 

21 

16 

0.2 

27 

15 

13 

0.5 

19 

6.5 

9 

1.0 

13.0 

0 

6 


Recall that the initial loop SNR is set equal to 1/^ FFT , the inverse of the variance of the FFT- 
estimated phase at the time the loop is closed. For the acquisition time plots presented in this section, it 
is assumed that S N Rfft is 12 dB, and that the tracking loop is closed at t c i t ioop = 'Rfft- In this case, 
the initial SNR L 's are listed in Table 4. For the SSL, 1/<^ FFT is less than the minimum loop SNR of 
10 dB required to avoid cycle slipping; therefore, the loop will be started with a bandwidth such that the 
initial loop SNR is 10 dB, as described in Section IV.F. By comparing Tables 3 and 4, it is seen that the 
carrier loop starts out very close to or above the SNRLjin required. Hence, the carrier loop acquisition 
time is essentially 0, leaving only the FFT acquisition time. As larger degradations to the output symbol 
SNR are tolerated, the required SNR L j in are smaller and the loop acquisition times are shorter. In fact, 
for a degradation of 0.5 dB, the loop acquisition times for both the symbol and the suppressed carrier 
are 0 sec, leaving only the FFT acquisition time. For a 1-dB degradation, the loop acquisition time is 0 
for all three loops. 

The loop acquisition time, Tyjoop, is calculated as the time to reach the required final bandwidth, 
Ri./iri) corresponding to the required SNRljuu Eq. (41), plus 1/2 BLjin sec. The latter term is an 
approximation of the additional time the loop transients will require to die out after the final window 
narrowing by an extra factor of two at the end of the narrowing process (Section IV. E). Hence, 
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1 


1 


1 


(46) 


Tloov — 


2 (dr/dt) \B lJin B L}init J " 2 B LJin 


where 


BLJin 


' 1 1 1 
c w c \ \JSN R.LJin 2?r 
< 1 

, <$SNR LJin { 2 tt )2 


SCL and SSL 
suppressed carrier loop 


and 


°J>FF 


BL t init — \ 


27T 


4>fft 


SCL and SSL 


suppressed carrier loop 


Table 4. Initial loop SNRs when SNRfft = 12 dB. 


Loop 

SNE-L^nit — FFT > ^B 

SCL 

14.8 

SSL (symbol phase 

8 (10 dB used to avoid 

estimation by interpolation) 

cycle slips) 

SSL (FFT frequency and 

8.8 (10 dB used to 

phase estimation) 

avoid cycle slips) 

Suppressed carrier loop 

14.8 


Acquisition times vary with the types of loops: one-arm versus two-arm and different loop orders. 10 
In this section, T sc iaC q, T$ym,acq-> a nd T carr > acq are presented for the high carrier phase-noise case. In 
Figs. 7(a), 8(a), and 9(a), T SC}(lcq is plotted for final degradations of 0.1, 0.2, and 0.5 dB, respectively. 
Similarly, in Figs. 7(b), 8(b), and 9(b), T sym ^ acq is plotted. The suppressed carrier acquisition time plotted 
in Fig. 7(c) is presented as T carr ac? = T carr>FFT (with the assumption that T carr j oop & 0 sec) for final 
degradations of 0.1, 0.2, and 0.5 dB. 

In these figures, Ty )acq is a function of Pd/N 0 and E s /N 0 (or i? 5ym ). Hence, families of curves of 
versus Pd/N 0 are plotted for various values of E s /N 0 . Note the decrease in the acquisition times 
for larger final degradations. 

The total acquisition time of the algorithm, T acq) for the final degradation of 0.1, 0.2, or 0.5 dB from 
each loop, is estimated as in Eq. (41). For example, for P F /N 0 = 14.6 dB-Hz and R S ym = 16 symbols per 
sec ( E s /N 0 = 2.56 dB) expected when the Galileo spacecraft is at JOI, the acquisition time, T acqiJO u for 
a final degradation of 0.2 dB is 


Tacqjoi = A/AA{26.4, 33.1} + 5.9 = 38 sec 


10 A complete set of acquisition times for the different loops is in C. Buu, S. Stephens, and M. Aung to J. Berner, "Acquisition 
Times of the Fast Acquisition Algorithm,” JPL Interoffice Memorandum 3338-94-052 (internal document), Jet Propulsion 
Laboratory, Pasadena, California, June 4, 1994. 
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- PREDICTED 
ACQUISITION 
TIME 

o SIMULATED 
ACQUISITION 


Fig. 7. Acquisition time to reach 0.1-dB degradation in output symbol SNR due 
to phase jitter in each loop: (a) subcarrier T sc gc g, and (b) symbol 7" S ym,acq’ 
and for (c) suppressed carrier ^carr,acq^° reac ^ or 0*5”dB degradation. 





The simulated acquisition times are also included in Figs. 7-9, shown in discrete points. In the 
simulations, SN RLjin padded by 1 dB was used. Relatively good agreement is seen between the predicted 
and the simulated times. The seemingly larger divergence in the plots for a 0.5-dB degradation is largely 
due to the vertical scale of the plot. In Figs. 7(c) and 9(b), no simulation points are included as T acq — 
Tfft- 




Fig. 8. Acquisition time to reach 0.2-dB degradation in the output SNR due to phase 
jitter in each loop: (a) subcarrier T sc gcq and (b) symbol T gym acq . 


VIII. Conclusion 

The fast acquisition algorithm developed for the Galileo S-band mission that is expected to encounter 
low data rate, low SNR, and high carrier phase-noise conditions is presented in this article. Multiple 
schemes resulting in faster acquisition times are introduced. Implementation of the algorithm in the 
Block V Receiver is described. Simulation work supporting the algorithm is discussed. Finally, the 
acquisition times of the algorithm are presented. In the expected case for the Galileo spacecraft at 
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in ii linn 


JOI {Pd /N o = 14.6 dB-Hz and R, ym = 16 symbols per sec), the acquisition time to attain a 0.2-dB 
degradation from each loop to the output symbol SNR is estimated to be 38 sec. 



— , dB-Hz 
N o 


Fig. 9. Acquisition time to reach 0.5-dB degradation in the output symbol SNR due to 
phase jitter in each loop: (a) subcarrier T sc aC(J an d (b) symbol T S y m gC g . 
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Appendix 

Frequency Interpolation and Phase Estimation From the FFT 


In this appendix, the technique adopted for estimation of the frequency and phase of an exponential 
tone from its FFT is discussed. Expressing the input exponential tone signal, s(m), sampled every T s> fft 
sec as 


s 


(m) = AS 2 **"'***) | 


t — mT w> FFT 


(A-l) 


the technique for estimating the frequency and phase (/ tn and <f> 0) respectively) is described. 


A. Frequency Estimation Using Sine Interpolation 

An A^px-point complex FFT of the signal s(m) will have a peak at bin number k pea k . By interpolation 
between the complex FFT values of the peak bin and two bins about the peak, n bins away, the frequency 
fin can be estimated as }fft • The interpolation assumes an underlying sine function relation between 
the three points (rather than a simple quadratic one), resulting in an improved estimate of frequency. 
The interpolated FFT estimate, f f ft > is 


f FFT = 


kpeak T A& 1 

Nfft T S) fft 


where 

A k = frequency bin interpolation 
kp ea k = bin number of the peak tone 
Nfft — FFT size 
T S] fft — FFT sampling time 

The Ak is calculated as 
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where 


Ak = n x Re 


peak mod N j? ft ^{Jtpeak ^)rr 


W FFT 


y Z(kp Ca k+n) 


modN p prp 


^(kpeak-n)modN FFT T (7 7)^ 


peak 


(A-2) 


Zk ~ complex FFT value of the Arth bin 

y — e j 2*nfa _ j 

NpFT 

Ot = — 

Ndata 

— zero-padding factor, may not be an integer 
n = MAX | Floor |],l} 

When a = 1, Eq. (A-2) becomes in the limit 


A k = Re 


Z Z 

“(hpc .fc + 1) mod W p prp ~ ^(*,..1-1) m od Wpp'j' 


Zn- 


(£j>ea*"l"l)rr»o<i;Vjrj; 


+ 3 


(^P 


_ 1 ) m od 


-2Z k 


pea k 


(A-3) 


B. Phase Estimation 

The phase 4> 0 of the signal at t = 0 can be estimated from its FFT as 4>fft ■ 


where 


&FFT = A <^FFT “ 2kJfFT S ’^ FT 


A<f) F FT — Arg 


[2* 


— 27rAA 


(■ Ndata - 1 ) 
2<a Ndata 


(A-4) 


(A-5) 


and T s fft is the sample rate of the FFT input signal. 
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